Atomic population distribution in the output ports of cold-atom 
interferometers with optical splitting and recombination. 

E. O. Ilo-Okeke and Alex A. Zozulya 

Department of Physics, Worcester Polytechnic Institute, 
100 Institute Road, Worcester, Massachusetts 01609, US^ 
(Dated: June 21, 2010) 

Abstract 

Cold-atom interferometers with optical splitting and recombination use off-resonant laser beams 
to split a cloud of Bose-Einstein condensate (BEC) into two clouds that travel along different 
paths and are then recombined again using optical beams. After the recombination, the BEC in 
general populates both the cloud at rest and the moving clouds. Measuring relative number of 
atoms in each of these clouds yields information about the relative phase shift accumulated by 
the atoms in the two moving clouds during the interferometric cycle. We derive the expression 
for the probability of finding any given number of atoms in each of the clouds, discuss features 
of the probability density distribution, analyze its dependence on the relative accumulated phase 
shift as a function of the strength of the interatomic interactions, and compare our results with 
experiment. 
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I. INTRODUCTION 



Using wave-like properties of atoms for atomic interferometry has been a subject of intense 
and extensive study [l| . Atoms are sensitive to electromagnetic fields due to their electric and 
magnetic moments, their mass allows them to be deflected in the gravitational field thereby 
making them attractive in the measurements of inertial forces 2|. Atom interferometers have 
been used to measure gravitational constant jsl, acceleration 4-6|, electric polarizability [3] 
and fine structure constant [8] to very high accuracy. 

The technical realization of neutral-atom interferometers took some time as compared 
to their electron- and neutron-based counterparts. Part of the reason is that atoms have 
large mass, resulting in a smaller de Broglie wavelengths for the same velocity. Also, neutral 
atoms cannot easily propagate in dense matter unlike, e. g., neutrons, and therefore require 
more ingenious ways to coherently split and diffract the atomic beam. The first atom 
interferometer {q] realized in a double-slit diffraction geometry, worked with a stream of 
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ll| used 



supersonic gaseous atoms and used mechanical gratings [9|. Later experiments 
standing light wave to coherently diffract the atomic beam. The standing light wave is 
formed using a laser beam that is detuned from atomic resonance to avoid spontaneous 
emission, and is retro-refiected by a mirror. The spatially-varying envelope of the standing 
wave creates an effective optical potential acting as a diffraction grating for atoms that can 
be used to split and recombine an atomic beam. 



Another technique j^, ll2|] for diffracting an atomic beam exploits Raman transitions 



between two hyperfine ground states of an atom via a third quasi-excited state. The laser 
pulses (often called Raman pulses), consist of two counterpropagating light beams with 
frequencies which are different by the Bohr transition frequency between the two hyperfine 
states. Absorption of a photon from one laser beam and stimulated re-emission into another 
one in this case is accompanied by a transition between the two hyperfine states. 



The use of Bose-Einstein Condensates (BECs) 13] in atom interferometers is appealing 
for many reasons. BEC has narrow momentum distribution that minimizes the spread in 
momentum during the splitting and recombination of the atomic cloud and reduces the 
expansion of the condensate during propagation. Ultracold BEC can be easily manipulated 
and confined in a very small area on an atom chip Finally, BEC has large coherence 
length allowing for good fringe contrast and helping to determine any offset phases more 
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accurately. Since the first experimental demonstration of interference between two different 
Bose condensates 15|, several experimental techniques for the manipulation of BECs and 
different BEC-based interferometric geometries have been proposed and demonstrated 



16|-ll9|. 



In trapped-atom interferometer geometries [16|, llTI] the BEC is kept in a trap confining 
the atomic cloud in all three dimensions. This trap is dynamically split into two double-well 
traps using a standing light wave, to create two arms which were physically separated in 
space. After some time, the trap is switched off allowing the condensates in each arm to 
fall, expand and interfere. 

In guided-wave interferometers the BEC is kept in a waveguide. The condensate is tightly 
confined in two transverse dimensions but allowed to propagate along the third dimension 
. The waveguide potential along this guiding dimension is typically weakly parabolic 
either because of difficulty of completely canceling magnetic field gradients or by design. 
A typ,ca. example is the M,che.so„-type smgle reflect,™ atom interferometer reaped m 
Ref. p^. In this interferometer, the BEC cloud tpo is initially at rest in a waveguide. 
Splitting pulses consisting of a pair of counterpropagating laser beams detuned from atomic 
resonance and acting as a diffraction grating, are incident on the cloud. These pulses split 
the condensate into two harmonics, and ip-y moving with the initial velocities ±fo, 
respectively. In a single reflection interferometer, the directions of propagation of these 
harmonics are reversed at time T/2 (where T is the duration of the interferometric cycle), i.e., 
in the middle of the cycle with the help of a reflection pulse. The harmonics are then allowed 
to propagate back and are recombined when they overlap again using the same optical pulses 
that were used to split the original BEC cloud. After the recombination, the condensate is 
in general in a superposition of ipo, tp+ and ip- with the relative amplitudes depending on 
the amount of the accumulated phase shift between the arms of the interferometer acquired 
during the cycle. 



In a double reflection interferometer 
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2l|, the optical reflection pulse is applied twice 



at times T/4 and 3T/4. After the first refiection pulse, the harmonics change their direction 
of propagation and start moving back. They pass through each other, and exchange their 
positions by the time 3T/4. The harmonic that was on the right at T/2 is now on the left 
and vice versa. The second refiection pulse applied at 3T/4 again reverses the directions of 
propagation of the harmonics and, finally, they are recombined at time T. 
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Authors of Refs. 
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23| investigated interferometric geometry that does not rely not 



on reflecting optical pulses but instead uses gradient of the confining waveguide potential 
for reversing direction of propagation of the BEC harmonics. In this "free oscillation" 
interferometer the moving BEC clouds propagate in a parabolic confining potential. They 
slow down as they climb the potential, stop at the their classical turning points after one 
quarter of the trap period (T/4) has elapsed, and turn back. At T/2 the clouds meet at 
the bottom of the potential, reach again their turning points at 3T/4 and are recombined 
at time T. The duration of the interferometric cycle is thus equal to the oscillation period 
of the parabolic longitudinal waveguide potential T. 

A Mach-Zehnder-type atom interferometer is another waveguide interferometer that uses 



BEC 
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22| . This interferometer shares the same operating principles with the above- 



discussed Michelson interferometer. The difference between the two interferometers lies in 
their splitting techniques. In the Mach-Zehnder-type atom interferometer the two coun- 
terpropagating waves used for a 7r/2 splitting pulse are frequency-shifted with respect to 
another resulting in a traveling optical potential. This tt /2 pulse transforms the original 
BEC at rest into two clouds of equal amplitude. One of these clouds remains at rest and 
the other propagates with velocity Vq. A tt pulse in the middle of the cycle stops the moving 
cloud and brings the one that was a rest into motion . Finally, the second tt/2 pulse applied 
at the end of the cycle recombines the two clouds. 

Both in trapped-atom and guided-wave interferometers, the interference fringes depend 
on the relative phase accumulated by the atomic clouds in different arms during the inter- 
ferometric cycle. Apart from the accumulated phase shift induced by fields or interactions 
of interest during the experiment, unwanted phase may be accumulated due to confinement 
effects and interatomic interactions resulting in a decrease in visibility of the interference 
fringes. For example, in the single-pass interferometer [l4|, during the propagation the outer 
edge of each cloud feels a higher potential than the inner edge (the outer edge is the leading 
edge in the first half of the cycle when the clouds move away from each other and the trail- 
ing edge in the second part when the clouds move toward each other). The outer edge thus 
accumulates a larger phase than the inner one. During the recombination, the outer edge 
of one cloud interferes with the inner edge of another and the phase difference accumulated 
due to the presence of the confinement potential leads to a coordinate-dependent residual 
phase across the clouds after the recombination. Another mechanism for phase accumula- 
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tion is due to mutual interaction of two BEC clouds ip- and ip^ when they spatially overlap. 
During the separation, the inner edge of one cloud interacts with atoms in the other cloud 
until it has traversed the entire length of the other cloud, while the outer edge of each cloud 
hardly interacts with any atoms in the other cloud (and similarly during the recombination). 
As a result, the inner edge accumulates larger phase than the outer edge. The above two 
contributions have opposite signs but different magnitudes so the net coordinate-dependent 
phase is not zero. Still another mechanism of accumulation of the unwanted phase is due 
to the fact that the velocities ±v of the moving harmonics 1^+ and ip- during the reflection 
are different from their initial velocities ±fo due to the influence of the confining potential 
and the interatomic interactions. As as result, the reflection pulses (which are formed by 
the same pair of counterpropagating beams as the splitting/recombination pulses) are not 
exactly on resonance and do not exactly reverse the clouds' velocities (the direction of prop- 
agation of each of the clouds does change but the speed before and the after the reflections 
is different). 

Theoretical analysis of the sin gle- and double-reflection interferometer geometries has 



been carried out in Refs. 
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24 



26l |. According to the analysis of Refs. 2l|, |26|], sym- 



metric motion of the two clouds in the double-reflection geometry partially cancels the 
velocity errors imposed by the reflection pulses and the phase imposed by the confining 



potential 
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26j. This conclusion has been confirmed experimentally in Ref. 2l|]. The 



"free oscillation" interferometer provides even greater degree of cancelation of the unwanted 
coordinate-dependent phase since it does not rely on reflection pulses and does not suffer 
from the velocity mismatch effects. Recent experiments 2l|, |22] where the atomic clouds 
were allowed to be reflected from their classical turning points instead of using reflection 
pulses to truncate their motion, confirmed more accurate cancelation of the unwanted phase. 

An additional mechanism that could lead to dephasing of the interference fringes is due 
to atom-atom interactions within each of the two clouds. These interactions result in the so- 
called phase diffusion in the BEC clouds 27H3l| . The aim of the present paper is to analyze 
effects of the atom-atom interactions within each of the BEC clouds on the operation of cold- 
atom interferometers using optical pulses for splitting and recombination of the condensate. 
The basis physics of the interferometric cycle in the presence of the interatomic interactions 
can be described as follows: The BEC which is initially in a number state with atoms is 
split by the optical pulses into two clouds moving in opposite directions. After the splitting 
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the system is in an entangled state with each cloud being a linear superposition of number 
states. This superposition is concentrated around the mean N/2 with the relative uncertainty 
of the order of 1/ V^. In the presence of interactomic interactions, each number state evolves 
with different rate resulting in the accumulation of relative time-dependent phases between 
the different number states. The recombination process is sensitive to these phases and thus 
the interactions should influence the contrast of the interferometric fringes. Specifically, we 
shall derive the expression for the probability density of observing any given number of atoms 
in each of the three output ports of the optical beamsplitter/recombiner and analyze it in 
different limits. Both the splitting and the recombination (detection) of the BEC are treated 
rigorously due to a large difference between the characteristic momenta of the condensate 
and the momentum imparted to the atoms by the optical beams. Thus our analysis does 
not rely on any ad hoc assumptions about the process of detection and can be considered an 
analysis from the first principles. The previously-discussed effects resulting in a coordinate- 
dependent phase due to confinement and the interaction between the overlapping clouds are 
not accounted for by the present analysis, i.e., we assume that most of the time the clouds 
are spatially separated. 

The remainder of the paper is organized as follows. In Sec. [Ill we discuss operation 
of the optical beamsplitter/recombiner and obtain the expression for the state vector of 
the system after recombination. In Sec. IIIIl we derive the expression for the probability 
density. Characteristic features of the probability density including the mean, the standard 
deviation etc., are calculated and the discussed in Sec. |Vl The results are compared with 
the experiment in Sec. IVII 

II. OPTICAL BEAMSPLITTER OPERATION 

Consider BEC cloud ipo at rest in a confining potential before the beginning of the inter- 
ferometric cycle. As was discussed in the introduction, the splitting optical pulses transform 
the initial BEC cloud into two clouds ip±i moving in the opposite directions. The clouds 
are allowed to evolve during the time T and are subject to the action of the recombination 
optical pulses (which are identical to splitting optical pulses). After the recombination, the 
atoms in general populate all three clouds tpo and iIj±. The relative population of the clouds 
depends on the phase difference between the clouds il)± acquired during the interferometric 
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cycle. 

Let 6q, b^_i and bl_i be operators which, acting on a vacuum state, create an atom be- 
longing to a cloud at rest and moving to the left or right, respectively. 

The many-body Hamiltonian describing atomic BEC during the interferometric cycle is 
of the form 

H = J dr¥{r,t)Ho^{r,t) + Y j dr¥{v,t)^^{r,t)^{r,0{r,t), (1) 

where Hq is a single-particle Hamiltonian and Uq = Airh'^asc/M, with M being the atomic 
mass and a^c the s-wave scattering length. The single-particle Hamiltonian Hq accounts for 
the confining potential for the BEC and also includes effects of the environment resulting in 
different rates of evolution for the BEC clouds propagating in opposite directions. Finally, 
the quantity ^(r,t) is the field operator 

^(r,t) = 6+iV^i(r,t) + 6_i^_i(r,t), (2) 

where ip±i{r,t) are wave functions of the BEC clouds moving to the right and left, re- 
spectively. The functions ip^^ cire shifted versions of each other propagating in opposite 
directions. 

Using Eq. (E]) in ([T]) results in the following Hamiltonian describing effects of the envi- 
ronment and the interactomic interactions: 

W 

^e// = -y (^+1-^-1) +^(^+1 + ^-1) • (3) 

Here n±i = bl_^b±i are the number operators, W is the relative environment-introduced 
energy shift between the right- and left-propagating clouds, and 

g = ^ j dr\^^,\^ = ^ j dr\^_,\^ (4) 

is the coefficient characterizing strength of the interatomic interaction. The Hamiltonian ([3]) 
neglects effects due to overlap of the right- and left-propagating clouds assuming that most 
of the time the clouds are spatially separated. 

The state vector of the system at the beginning of the interferometric cycle before the 
splitting pulses is given by the relation 
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where |0) is the vacuum state and is the total number of atoms in the BEC. 

Sphtting and recombination pulses couple the operators b±i and 60 according to the rules 



(see Appendix in Ref. 



251 ] for the corresponding splitting/recombination matrices): 



1 e"'^ 1 

bo (b+i + 

1 e"*^ 1 

where /3 is a phase factor. 

A single-atom state is transformed by the splitting pulses as 



(6) 



bl\0)^^{bl, + bU)\0), 



so the product state vector of the A^— particle system Eq. ([5]) after the splitting acquires 
the form 

N y^r\ 1/2 



where 



|n+,n_) 



■|o) 



(7) 



(8) 



is the state with n+ atoms traveling to the right and n_ to the left, respectively, and 
(^) = A^!/n!(A^ — n)\ is the binomial coefficient. 

Time evolution of the state vector is governed by the Hamiltonian (j3]): 

ft 



= l"^ split) exp 



{t/h) I H^ffdt' 

Jo 



States |n+, ?T,_) given by Eq. ([8]) are eigenstates of the Hamiltonian (|3]) with the eigenvalues 



W 

E{n+, n_) = — -{n+ - n_) + g{n^ + n_). 



(9) 



The state vector of the system at time T immediately before the recombination is thus given 
by the relation 



|^(T)) = ^^> I le 



^\b\,\ [bl. 



N-n 



|0). 



(10) 



Here 



1 



T 



e = - I dtw (11) 



is the environment-introduced accumulated phase difference between the right and the left 
clouds and 



« = i£*9 (12) 

is the accumulated nonlinear phase per atom due to interatomic interactions. The recom- 
bination pulses act on |\1/(T)) in accordance with Eq. ([6]). The resulting state vector of the 
system after the recombination has the form: 

N 



n=0 

N-n 



where we have omitted irrelevant phase term exp(— 

III. PROBABILITY DENSITY 

The state with atoms being in the cloud moving to the right, n_ in the cloud moving 
to left and hq = N — — n_ in the cloud at rest, is described by the state vector 

{^ur ("-0"" w"° 



The probability of this outcome after the recombination is given by the modulus squared of 
the probability amplitude (n+, n_, no|\l/rec)- Using Eq. ( fT3l) . this probability amplitude can 



be written as 

N 



(„| (to)"" (6-1)"- (6+1)"* 



^ie{n~~N l'i)-2ii(n-N I2f 



\fno\ ^yn^l 
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Equation (fT5|) can be recast as 



[n+,n_,no\ 



N 

N/2)-2ii(n-N/2f 



n=0 

1 (9"o 9"- (9"+ 



n+ 

i/3 



^4 r^+^^n 10). (16) 



2 72 7 \^ 2 V2 

The product of two terms in parentheses can be represented as the double sum 



N-n 



(17) 



1=0 j=0 

The derivatives with respect to Bq in Eq. f|T6l) will select only the term with i + j = fiQ from 
this sum, yielding 



n=0 

where 



fN\fN\-^ n N-n / \ /at \ 

\n)\no) (n\(N-n^ 

i=0 j=0 

and 

a = /3no + (7r/2)(n_-n+). (20) 

The unwieldy expression ( fT9l) can be written in a much more manageable form by evaluating 
Eq. ( JT6l) for = 0. In this case, summation in Eq. (JT6l) can be readily carried out and, after 
differentiation, Eq. (fT6l) results in the expression 



K,„_,n„|*.„) = e-y-^-^(^cos-j (^sm-j . (21) 

Comparison of Eq. f|T8|) and Eq. ( 12T|) shows that 

5:e-^/(n,no) = e^^^/^(cos^) (sin^ (22) 

n=0 ^ / V / 
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immediately yielding 

/(n, no) = - r dxe'"(^-2") (cos (sin xf""" , (23) 
Jo 

Using Eq. ( !T8|) . we can write the probability density 

P(n+,n_,rao) = | (n+, n_, rao|^rec)P 
as the product of two functions: 

P(n+, no) = P±(n+, n_, no)Po(^o, 6^, 0, (24) 

where 



^ (iV-no)! 
^ n+!n_!2^-"o 

and 



(25) 



Po(^o,^,0= (^J |S(no,^^,Ol'- (26) 
with the function S(no0,^) given by the relation 

N 
n=0 

The function P± describes the probability of observing n+ and n_ atoms in the right and 
left clouds, respectively, for any given number = N — (n+ + n^) atoms in the central 
cloud. This function is independent both on 6 and the nonlinearity ^ and is normalized to 
one: 



N—no N—no 

n^\{N — no — n+)! 



V P - V (^-^°)'^"°"'' - 1 (OR) 



n+=0 n+=0 

With the use of Stirling's approximation in Eq. ( l25l) . P± can be simplified to: 



2 

exp 



vr(iV-no) 



|t^-,„,-H„_^iV-,,. (29) 



The function Po(no, 6,^) describes the probability of observing no atoms in the central cloud. 
The effects of both the external phase 9 and the nonlinearity ^ are contained in this function. 
It is also normalized to one: 

N 

J^Po(no,^,0 = l- (30) 

"0=0 
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The function S ( 127|) satisfies the symmetry relations 

E(no,-^^) = (-l)^-"«S(no,^), 

E{no,7r -6) = ^{N -no,e). (31) 

The probabihty function Pq given by Eq. fl26|) is periodic in ^ with the period 27r. Relations 
( l3T]) alow us in the following to restrict our analysis to the values of 6 lying in the interval 
0<9<n/2 since 

Po{no,7r-e) = Po{N-no,e). (32) 



IV. EVALUATING THE PROBABILITY DENSITY FUNCTION R 







The exact expression for the probability density distribution function Pq given by the 
equations f l2B]) . fl2^ and fl27|) . is relatively complex and does not lend itself readily to an 
easy interpretation. In the following we shall transform and simplify it to make it more 
amenable for the subsequent analysis. 

The integral I{n, hq) given by Eq. fl23|) can be evaluated in the complex plane by the 

method of steepest-descent to yield: 

1 



I{n,no) 



exp 



fAT M /I ^0 , 1 {n-N/2f 
A/ — nn m 1 / 1 h nn m 1 / 



X 



AN- 



-2n) arccos ^uq/N _j_ ^ j^^A''—nog—i{A''—2n) arccos 



/TV 



(33) 



Using Eq. (l33l) in the expression for S (1271) . approximating summation by integration and 

evaluating the integral, we get: 

1 



S(no) 



exp 



X 



VI + 2iiN 



(AT -no) In ^/l - ^ + "-0 I^a/^ 



where 



f arccos y/no/N =f 6'/2 



(34) 



(35) 



1 + 2i^N 

Finally, using the Stirling's approximation and Eq. ( |34l) in Eq. (l26ll results in the expression 
for the probability density Pq: 

1 



Poino) 



N 



V27r(l + 4^2;V2) y no{N - uq) 



(36) 
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Equation (136|) is not applicable at the two end points, no and uq = N, where it has to be 
replaced by the expressions 



Po(0) = , exp 
Po{N) = , ^ exp 



-2(71 -9fN 
1 + 



V. FEATURES OF THE PROBABILITY DENSITY 



(37) 



Expressions ( l25i) and ( 136|) for the probability density functions P± and Pq give the prob- 
ability P(?2+,n_,no) = P±Po of observing any given number of atoms in the three output 
ports (three atomic clouds) of an optical beamsplitter. The function P± describes the prob- 
ability of observing n+ and atoms in the right and left clouds, respectively, for a fixed 
number hq = N — (n+ + n_) atoms in the central cloud. This probability is the Gaussian 
distribution f l29|) with the average values of n„ and n+ given by 

(n_) = (n+) = i(iV-no) (38) 

and the standard deviations 

An_ = An+ = ^A/iV - riQ. (39) 
The numbers of atoms in the right and left clouds are anti-correlated: 

Cov(n+, n_) = {n+n_.) - (n+) (n_) = -^(iV - ^o). (40) 

The probability to find nQ atoms in the central cloud is given by the function Po(no,6',,^) 
Eq. (l36l) . The dependence of this function on its arguments is not trivial, so we start the 
analysis by evaluating the expectation value of the atoms on the central cloud (uq) and the 
standard deviation Auq. 

The function Pq is proportional to the modulus squared of the sum of two terms: Pq oc 
le-''- + (-l)^-"»e-''+ , where rjzp are given by Eq. (l35l) . The relative phase difference 
between them, as a function of uq, changes rapidly due to the multiplier (—1)"''. Thus, 
the interference terms can be neglected in calculating both the mean and the standard 
derivation: 

N 



(Tin) ~ f dVn ! fp-SRe,,. , -2Re^+\ 
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Evaluation of the above integral yields 

(no) = y [1 + exp {-2eN) cos 6] . (41) 

Similarly, the standard deviation evaluates to 

(Ano)^ = ^ [1 - exp (-4^^iV)] [l - exp (-4^^N) cos 2^1 . (42) 
8 

To understand these results, we shall look at the dependence of the function Pq on the 
number of atoms uq for different values of the strength of the interatomic interactions ^. At 
relatively small values of ^ such that ^ <^ 1/\/N, the term exp(— ?7_) in the expression (j36l) 
for the probability dominates the second one. The probability Pq is then a simple Gaussian 

21 



Po oc exp 



2N (arccos ^n^/N -6/2 



1 + 4^2 JV2 



(43) 



with the maximum located at = N cos"^ 6/2. This situation is shown in Fig. [T] The 
two curves in this figure are plots of the function PqIuq) given by Eq. fl36l) versus no for two 
difference values of the interatomic interaction strength C,- Both curves correspond to the 
same value of the angle 6. The most noticeable feature of Fig. [T]is the increase in the width 
of the probability distribution with ^. This behavior is explained by Eq. (j4T|) . which in the 
limit ^ ^ 1/\/N reduces to: 



N 



Ano ^ ^ sin ^ ^/TTWn^- (44) 

For very small values of (^ ^ l/^)? the influence of the interatomic interactions on the 
operation of the beamsplitter is negligible. The relative standard deviation of the number of 
atoms in the central cloud is inversely proportional to the square root of the total number 
of atoms in the system: Auq/N oc 1/y/N. For ^ ^ <^ the width of the 

distribution linearly grows with the increase in ^. 

The mean value of for ^ ^ 1/ \/N reasonably closely corresponds to the position of 
the peak. Equation ( HT|) for (uq) in this hmit yields: 

(no) ^ -(1 + cos^). (45) 

As is seen, no depends on 6 but not on ^. 
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For larger values of ^ ^ 1\/N, the width of the probability distribution function Pq 
becomes of the order of the total number of atoms in the system N. The two terms exp(— 
and exp(— ?7+) in Eq. f l36|) are now comparable in magnitude. The transition to this limit is 
shown by Fig. |2] and Fig. |3l Black regions not resolved in Figs. |2] and [3] correspond to rapid 
spatial oscillations with the period 2. These oscillations are clearly seen in Fig. HJ which 
shows part of Fig. |3]for a narrow range of values of riQ. The oscillations are caused by the 
interference between the two terms in Eq. ( 136|) . As the magnitude of ^ approaches 1/\/N, 
these terms become comparable in magnitude. Because of the nearly vr— phase change 
between the two terms every time uq changes by one due to the factor (—1)"° , the two 
terms consecutively add either in phase or out of phase when one steps through different 
values of no. Along with rapid spatial oscillations, both Fig. [2] and Fig. [3] demonstrate 
oscillations of the envelopes at a much slower spatial rate which are more pronounced for 
larger values of the interaction strength. These oscillations are due to the fact that the 
relative phase of the terms exp(— r^^) and exp(— r^^) in Eq. (!36|) changes with uq. The nodes 
in Fig. [3] correspond to the value of this relative phase being equal to or a vr and an 
antinodes have the phase shifted by ±7r/2. 

Figs. [2] and [3] indicate that the probability Pq and, as a consequence, (uq) and Auq, 
become less sensitive to changes in the environment-introduced angle 6. This fact is graph- 
ically illustrated by Figs. |5] and [6] showing the average value of the number of atoms in the 
central cloud (rio) and the standard deviation Auq versus 6 as given by Eqs. fj4T]) and fj42l) . 
respectively. Fig. [5] demonstrates that increased interactomic interactions eventually lead 
to the loss of contrast of interference fringes. Additionally, larger interactomic interactions 
cause larger shot-to-shot fluctuations in the number of atoms in each of the three output 
ports, as is seen from Fig. O 

The loss of contrast of the interference fringes can be quantified by writing Eq. fj^Tj) as 

N 

{no) = -{l + Vcose), (46) 

where 

V = exp {-2^^N) (47) 

is the fringe contrast. Figure [7] shows the fringe contrast V Eq. (H71) as a function of ^ and 
demonstrates that the values of ^ approaching 1 / \/N result in a washout of the interference 
fringes. 
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VI. DISCUSSION 



The value of the accumulated nonlinear phase per atom due to interatomic interactions £ 



Eg. ([T21) depends on the volume of the BEC clouds (cf. Eq. (jl])). Experiments 
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Q,[l8|, 



20 
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22| to be discussed below, were conducted in parabolic traps with confining potentials 



of the form 

M 
~2 



Density profiles of the moving clouds are well described by the Thomas- Fermi approximation 

(recall that V'i are normalized to one). 

Immediately after the splitting pulses, the density profiles of the moving clouds are the 
same as that of the initial BEC cloud containing atoms and being in equilibrium in the 
confining potential Eq. ( |48|) . After the splitting, each moving cloud contains A^/2 atoms. 
The repulsive nonlinearity is no more balanced by the confining potential and the radii of 
both clouds start to oscillate. The maximum size of the oscillating clouds is the equilib- 
rium size corresponding to N atoms and the minimum size lies below the equilibrium size 



corresponding to A^/2 atoms. For estimates, we can take in Ea 



equilibrium size of a cloud with A^/2 atoms: = 2fi/Muf, where [321. l33| 



to be given by 



15 V"/..U„^"'' 



Uq = Anh'^asc/M, uj = {ujxOJyUizY^'^ and a = ^Jh/Muj. 

Evaluation of the constant g Eq. (jlj) yields g = {15Uo)/ {28'7iRxRyRz)- The accumulated 
nonlinear phase per atom due to interatomic interactions ^ Eq. ( IT2l) is then given by the 
expression 

(^'-{mV2y'\^y"^TN-:'", (51) 

where T is the duration of the interferometric cycle. 

The relative importance of interatomic interaction effects on the operation of the inter- 
ferometer is determined by the parameter P = ^\/N : 

P = Om(^^Y\loT)N~'/''>. (52) 
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Figure ([7]) shows that the contrast of the interference fringes decreases with the increase 
in P. The condition of good contrast can be somewhat arbitrarily stated as P < 1/2 (for 
P = 0.5, V = 0.6). 

Equation fl32]) shows that P oc Tu^^^N^^^^^ . The dependence of P on the total number 
of atoms in the BEC clouds is very weak, and so this parameter is primarily dependent on 
the duration of the interferometric cycle and the averaged frequency of the trap. 

Experiments by Wang et al [l^ were conducted using the Michelson geometry. The BEC 
consisted of about 10^ ^''Rb atoms [34]. The transverse and longitudinal frequencies of the 
trap were 177 Hz and 5 Hz respectively. The propagation time T was up to 10 ms. For these 



parameters and the value of the scattering length = 5.2 x 10^^ m 35|, Eq. (132]) yields 
P ^ 1.6 X 10~^. Thus, the interatomic interactions were not limiting the visibility of the 
interference fringes in these experiments. 



A similar experiment was performed by Garcia et al 20|, |2l| also in the geometry of a 



Michelson interferometer. In Ref. 



20|, a BEC cloud of about 10 Rb atoms has been 



produced in a trap with the frequencies of 6.0 Hz, 1.2 Hz and 3.0 Hz, respectively. The 
interferometric time T was about 40ms. Using Eq. ( l52l) . we can evaluate the value of the 
parameter P in the experiment as P ~ 10^^, which was too small to result in observed 
degradation of the contrast. The loss of visibility in the experiments |20| was attributed by 
the authors to spatial noise on the splitting beams and asymmetric splitting of the cloud 
due to the condensate's residual motion when it was loaded into the trap. At longer tines, 
the loss of coherence might have been caused by various noise sources. Similar results were 
reported in Ref. 21|, were the confinement frequencies were deliberately kept weak making 
the atomic density and this the interatomic interaction effects small. 



Horikoshi et al 
of atoms in Ref. 
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22| demonstrated a BEC Mach-Zehnder interferometer. The number 



22| was about 3x10^ and the radial frequency of the trap was fixed at 60 Hz. 
The experiments have been conducted for two different values of the axial frequencies and 
interrogation times T. At an axial frequency of Uz = x 17 Hz and the propagation time 
of the cloud about T = 60 ms, the parameter P = (f)^/N estimated using Eq. ( l52l) turns out 
to be about 0.38. For this value of P, Eq. P7|) gives the value of the fringe contrast about 
70%. The experimental value is 30% [22]. Similarly, for the axial frequency was 10.29 Hz 
and the interferometric time 97 ms, Eq. f l52|) gives the value of P ~ 0.5 corresponding to an 
estimated contrast of 58%. In this case no fringes were observed experimentally with a about 
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40% scatter of the data points. The authors of Ref. 22| conjecture that vibrations could be 
the main source of the loss of fringes in their experiments. The above estimates indicate that 
the interatomic interactions discussed in the present paper could be also partially responsible 
for the observed degradation of the interference fringe. 
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Appendix A: Mach-Zehnder-type interferometer 

In a Mach-Zehnder-type cold-atom interferometer, the optical splitting 7r/2 pulse trans- 
forms a BEG cloud at rest in a superposition of two clouds ipo and ip+i- The first cloud is 
at rest and the second one is moving. The clouds evolve during the time T/2 and are then 
subject to the action of a vr pulse. It stops the moving cloud and brings the one at rest into 
motion, i.e., transforms the ipo cloud into ip^i and vice versa. After additional evolution 
time T/2, the clouds are subject to a recombination tt/2 pulse. After the recombination, 
both ipQ and ip^i are in general populated. 

Analysis of a Mach-Zehnder-type interferometer parallels that given in the paper for 
the Michelson-type interferometer and is somewhat simpler because with the Mach-Zehnder 
interferometer there are only two output ports as opposed to three in the case of a Michelson- 
type interferometer. 

Let and 6^^^ be operators which create an atom belonging to a cloud at rest and moving 
to the right, respectively. The Hamiltonian is of the form (cf. Eq. ([3])): 

^e// = -y (^+1 -^o) +^ (^+1 + ^o) ' (Al) 

where the notations are the same as in Section [Tll 

The state vector of the system at the beginning of the interferometric cycle before the 
splitting pulse is given by Eq. 

The splitting/recombination tt/2 pulse couples the operators bi and bo according to the 
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rules: 

feo {bo + ih+i) , 

6+1 ^ i= {iho - 6+1) . (A2) 

For the tt pulse, similarly, one has: 

60 ib+i, b+i -ibo. (A3) 

Repeating steps of Section [Tll we arrive at the following expression for the state vector of 
the system after the recombination pulse: 

I ^ie{n-N/2)-2i^{n-N/2f 



X {bl-tbl,y {bl + zbl,Y'^^\0) (A4) 



The probability to have after the recombination hq atoms at rest and = N — no atoms 
moving is given by the modulus squared of the probability amphtude (no, N — nol'^rec)- Here 

Kn+)= ^/L ^ |0) (A5) 

V^o! Vn+! 

is the state with no atoms at rest and n,+ atoms moving, respectively. 

Repeating steps of Sec. Illlt the probability P{no) = |(no, — no\'^rec)\'^ can be written 

as 

N\ ,„ 2 



P{no,e,0= |S(no,^,OI , (A6) 

where S is given by Eq. (I27|l with the function /(n, no) in Eq. (1271) given by the relation 
(125]) . The probability P Eq. flA6p is thus completely identical to the probability Po Eq. fl2^ 
of Section IIIII All relations of Sections |IV] and |VT] equally apply to the case of the Mach- 
Zehnder-type interferometer. 
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FIG. 1: Probability function Pq versus no for ^ = and ^ = 3/N. For both curves, 9 = 7r/4 and 
N = 2000. Note that the abscissa axis range is from no = 1000 to 2000. 
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FIG. 2: Probabihty function Pq versus no for ^ = 0.2/y/N, 9 = tt/4 and AT = 2000. 
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FIG. 5: Normalized mean value of the number of atoms in the central cloud {no)/N versus 6 for 
= 2000. 
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FIG. 6: Normalized standard deviation Ano/N versus ^ for N = 2000. 
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FIG. 7: Interference fringes contrast F as a function of the strength of the interatomic interactions 
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